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Abstract. We improve the McLerran-Venugopalan model [1,2] by introducing a 
charge-density correlation function which is consistent with the observation that nu- 
cleons carry no net color charge. The infrared divergence in the transverse coordinates 
that was present in the McLerran-Venugoplan model is eliminated by the enforcement 
of color neutrality. 

The problem of extracting first principles predictions from the theory of strong 
interactions, quantum chromodynamics (QCD), is notoriously difficult, largely be- 
cause of the nonlinearity of the theory. The quanta of the gauge field, the gluons, 
themselves carry color charge and so serve as a source of additional gluons. Thus, 
any regime in which it is possible to actually compute some physical observable 
within the framework of QCD is of great interest. 

Such is the case with the McLerran-Venugopalan (MV) model [1,2]. What McLer- 
ran and Venugopalan realized is that for very large nuclei at very small values of 
the longitudinal momentum fraction, the number of color charges participating in 
the generation of the QCD vector potential is large. In this situation, the gluon 
number density may be approximated to lowest order by solving the classical Yang- 
Mills equations in the presence of the (classical) source generated by the valence 
quarks. To actually extract the gluon number density, we must average over the 
sources to obtain the two point correlation function for the vector potential. This 
approximation may be systematically improved by including the quantum correc- 
tions. The correlation function derived in Ref. [2] is highly infrared divergent at 
large transverse distances. In this talk, we will show that this difficulty may be 
ameliorated by forcing the nucleons to obey a color-neutrality condition [3]. 

Before proceeding with the main part of this talk, it is necessary to say a few 
words about our notation and conventions. We elect to work with light-cone coordi- 
nates, defined by x ± = (x° ±x 3 )/^/2. The components of a 4- vector will be written 
) , with vectors in the two-dimensional (transverse) subspace writ- 
ten with arrows. We choose a metric with the signature (— , +, +, +). Consequently, 
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cross sectional area of tube = d x 




FIGURE 1. A cartoon of a very large nu- 
cleus in its rest frame. In the labatory frame, 
Lorentz contraction causes all of the quarks 
and nucleons to pile up at essentially the same 
value of x~ . 
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FIGURE 2. Map of transverse momentum 
scales in the McLerran-Venugopalan model. 



we have x± = —x T , and a dot product which reads x • q = —x + q~ — x~q + + x • q. 
We will work in the light-cone gauge, A + = 0, where the intuitive parton model is 
realized [4]. 

Within the MV model [1,2], the Lorentz-contracted nucleus is treated as a pan- 
cake of color charge coming from the valence quarks. For a large nucleus, this 
means that a tube of cross-sectional area cPx will intercept a large number of va- 
lence quarks (see Fig. 1). If we restrict ourselves to small longitudinal momentum 
fractions, then all of the quarks in the tube effectively overlap. Because a large 
number of quarks contribute, the color charge is in a high-dimensional representa- 
tion of the group, and may be approximated by a classical charge density. 

To compute the gluon number density, we begin by solving the classical Yang- 
Mills equation for the vector potential as a function of the color charge density p. 
Next, the usual quantum mechanical average is approximated by an average over an 
ensemble of nuclear charge distributions to compute two-point correlation functions 
of the vector potential, (A°(x~, x)A b (x'~, x')). These correlation functions contain 
the gluon number operator; hence, they are be connected to the gluon number 
density. The result of the averaging process depends on the form chosen for the 
two-point charge density correlation function {p a (x~, x)p b (x'~, x')). Finally, we 
obtain the gluon number density by performing the appropriate Fourier transform 
of the position space correlator. 

Under what conditions is it legitimate to use the classical approximation we have 
just described? Schematically, the situation is as illustrated in Fig. 2. If the typical 
transverse momenta being considered is too small (below Aq CD ), then we are in 



the stong coupling regime, where we have no reason to expect classical methods to 
be valid. At larger transverse momentum scales, we have weak coupling. However, 
we cannot go to too large a transverse momentum, or else the requirement that 
we have a large number of quarks contributing in each patch d 2 x of nucleus no 
longer holds: the quantum graininess begins to become important. At the upper 
end of where the classical approximation is valid, the solution is QED-like, with 
a distribution function which is proportional to l/<f 2 . Somewhere in between is a 
cross-over scale where the non-Abelian terms in the equations of motion become 
important, and we require a solution which is summed to all orders. 

We now turn to an examination of the basic form of the classical solution, using 
the situation in electrodynamics as our guide. Consider a point charge e moving 
at the speed of light down the z axis. In this situation, the current is 

J+ = e 5(x')5 2 (x); J" = 0; J — 0. (1) 

Unlike in the (more familiar) Lorentz gauge where the /ith component of the current 
generates the /ith component of the vector potential, in light cone gauge only the 
transverse components of are non-zero: 

A'(x) = ^Q{x-)^. (2) 
This vector potential corresponds to the non-vanishing field tensor components 

These components of the field tensor correspond to transverse electric and magnetic 
fields of equal strength. There are no longitudinal fields in this limit. An observer 
sitting at some fixed position (6, x 3 ) would see no fields except at the instant when 
the charge made its closest approach, when a 5-function pulse would be seen. The 
magnitude of the pulse would be proportional to 1/b. The (unobservable) vector 
potential, which was zero before the passage of the charge, takes on a non-zero 
value for all times afterward. Since the fields vanish at these times, the late-time 
vector potential may be thought of as some particular gauge transformation of the 
vacuum. 

The overall features of the above description continue to hold when we switch to 
QCD, although the details are slightly altered by the presence of the non-Abelian 
terms in the field equations. The net effect of these terms is to color-rotate the 
source in a complicated fashion. Nevertheless, the chromoelectric and chromomag- 
netic fields are non-zero only at the instant of closest approach by the charge, and 
the vector potential switches from one gauge transform of the vacuum to a different 
one at this instant. 

We now outline the method used to compute the correlation function 
(A?(x~, x)A b j(x'~, x')). Full details of this part of the calculation are found in 



Ref. [2]. Essentially, what one does is to expand the vector potential in "powers" 
of the charge density p, and perform the averaging by doing all possible pairwise 
contractions using 



(p a (x-,x)p b (x'-,x')) 



5 ab p 2 (x~ 



5(x — x' )V(x — x) 



(4) 



In this expression p?(x~) is the color charge squared per unit area per unit thickness 
(i.e. per unit x~). We set up the calculation with a non-zero nuclear thickness to 
avoid ambiguities in the commutator terms of the Yang-Mills equations which would 
arise if we let p ~ S(x~) exactly. Refering back to Fig. 1, we see that the quarks 
in our "tube" of color charge typically come from different nucleons. Thus, we 
expect them to be uncorrelated, hence the dependence 5(x~ — x'~) in Eq. (4). The 
transverse dependence of the charge denisty correlator is given by T>(x — x'). In 
Ref. [2], it is argued that this should also be a delta- function, since we are restricted 
to length scales < Aqq D . However, as we shall see, doing so is not consistent with 
color neutral nucleons, and leads to severe infrared divergences in the correlation 
funtion. 

After doing all of the contractions, we resum the series to obtain the master 
formula 
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N c L(x - x') 



,N cX {x- 



-)L(x-x') 



(5) 



Eq. (5) depends on two new functions. The first, x( x ~ i x '~)i measures the amount 
of charge in those layers of the source which have already passed both of the points 
which we are comparing: 
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The appearance of this function may be understood by recalling that the value of 
the vector potential depends on whether or not the charge has yet reached its point 
of closest approach. Although the range of integration in (6) extends to x~ = — oo, 
in practice this is cut off by the form of p 2 (x~), which for a pancake-shaped charge 
distribution, should be non-zero only in a relatively small range near the value of 
x~ that corresponds to the position of the nucleus. 
The second new function appearing in (5) is given by 



L(x-x') = J d 2 ij d 2 i'V(i-i')[G(x-i)G(x' -f) 
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where 



This Green's function is the inverse of the operator V 2 (in two dimensions). The 
lack of a scale in our theory produces an infrared divergence in this function, which 
is signalled by the appearance of an arbitrary length scale A. Clearly, the infrared 
finiteness or lack thereof of L(x — x') is intimately related to the form chosen for 
— £')• O n the other hand, according to Eq. (7), L(x — x') vanishes when x = x', 
i.e. in the ultraviolet. Thus, at very short distances the nonlinear terms in (5) drop 
out and the behavior of the correlation function is the same as if we had considered 
a purely Abelian theory instead. 

If, as in Ref. [2], we take D(|*- £') = <5 2 (£ - we end up with 



L(x-x>)~^{5 2 (x-x')-5 2 (0) 



(9) 



Although the subraction term serves to remove the quadratic infrared singularity, 
a logarithmic divergence remains. Assuming that the arbitrary scale should be of 
order Aqcd on physical grounds, the authors of Ref. [2] obtain 
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This correlation function diverges like (x 2 ) x2 for large values of the separation. 
Of course, the bad behavior does not begin until the point x ~ Aq^ d , the point 
where we begin to mistrust our calculation anyhow. Unfortunately, because of this 
divergence, the Fourier transform of Eq. (10) does not exist for any value of q, real 
or imaginary. So it is difficult to see how to define the gluon number density more 
than qualitatively using this expression. 

The resolution of this problem lies in recognizing the importance of enforcing the 
observation that nucleons, when observed on a large enough length scale, should be 
color neutral. A consequence of the Gaussian averaging employed in the MV model 
is that the average color charge vanishes, (p a (x~,x)) = 0. However, we should also 
impose the (stronger) condition that 



J dx~ d 2 x p a {x~,x) = 



(11) 



for a nucleus-sized volume. If we integrate the charge density correlator (4) over 
all (x,x') and apply (11) we obtain a constraint on the transverse portion of the 
correlator V: 



J d 2 x V(x) = 0. 



;i2) 



Any correlation function which satisfies Eq. (12) is compatible with color neutral 
nucleons. Such a correlation function must contain an intrinsic scale, i.e. the mini- 
mum transverse length scale for which (12) becomes true. When we compute L(x), 



this scale will be imparted to the logarithms appearing in Eq. (7). This stongly 
suggests that the resulting correlation function determined from Eq. (5) will be 
infrared finite. Indeed, this is the case [3]. 

To illustrate the features of our improved treatment, we turn to a specific model 
of a large nucleus introduced by Kovchegov [5]. In this model, we view the nucleus 
as containing A nucleons of radius a distributed within a sphere of radius R. Each 
"nucleon" consists of a qq pair. By explicitly averaging over the allowed positions 
of the quarks, antiquarks, and nucleons we may explicitly compute the function 
T>(x — x'). We find that there are two types of terms. The first is generated 
when the position of two quarks (or two antiquarks) overlap. It is proportional to 
S 2 (x — x'), precisely the form for V employed in Ref. [2]. The second kind of term 
enters in with opposite sign and corresponds to the situation when a quark overlaps 
an antiquark. It is proportional to a smooth function of the separation. Since the 
focus of Ref. [5] was on very short distances, this term was neglected relative to 
the delta-function contribution. However, at somewhat longer distance scales, it is 
precisely this additional contribution which is required to satisfy Eq. (12). 

We have compiled a series of plots (Figs. 3-6) to aid in the comparison of our 
results using Kovchegov's model to the results of Ref. [2]. In addition to the uniform 
distribution of quarks, antiquarks, and nucleons employed by Kovchegov, we have 
also performed the averaging using Gaussian distributions. In preparing these plots, 
we have adjusted the nucleon size parameter a and (for the MV result) Aqcd so 
that the corresponding correlation functions (5) match in the ultraviolet limit. This 
requires ac = 0.464a;/ and Aqcd = 1.44a [ } 1 where ay is the nucleon size parameter 
for the uniform quark distribution. 

In Fig. 3 we have plotted the smooth part of the correlation function for the 
uniform and Gaussian cases, defined by writing T>[x — x') = 5 2 (x — x') — C (x — x') . 
In both cases, the bulk of the contribution comes from separations less than 2a. 2 
Therefore, in this and subsequent plots we have defined the dimensionless distance 
X = \x-x'\/{2a). 

Next, we present Fig. 4, which shows the trace of the correlation function (5) 
for all three cases. The Fourier transform of this quantity is proportional to the 
gluon number density. The Gaussian curve, not surprisingly, has a Gaussian tail 
for X > 2a, while the uniform curve vanishes identically in this region. On the 
other hand, the result from Ref. [2] runs off to — oo beyond Aq CD . 

We come to the gluon number density in Fig. 5. To define the MV curve, we 
have followed the suggestion made in Ref. [6] , and simply cut off the x integration 
at x — Aq CD . This is the source of the wiggles visible in Fig. 5. All three curves 
have the same overall shape, with a plateau at small values of q and a 1/q 2 fall off 
at large q. The number of zero momentum gluons clearly depends on how we have 
cut things off: the Gaussian model, which allows quarks to be (albeit with small 
probability) a large distance from the center of their nucleons has the most long 



2 ) In fact, C(x — x 1 ) vanishes identically for separations greater than 2a when a uniform distri- 
bution is assumed. 
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FIGURE 3. The smooth part of the 
two-point charge density correlation function 
in Kovchegov's model, defined by writing 
T>(x) = 5 2 (x) — C(x). The nucleon size param- 
eters have been chosen so that the resulting 
gluon number densities match in the ultravio- 
let limit. The two curves correspond to uni- 
form and Gaussian distributions of the quarks 
inside the nucleons. 
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FIGURE 4. The trace of the two-point vec- 
tor potential correlation function (5) in posi- 
tion space. The nucleon size parameters a and 
Aqcd have been chosen so that the resulting 
gluon number densities match for X — > 0. We 
have fixed the longitudinal coordinates at a 
place where N c \ = 50a -2 . Plotted are the re- 
sults of Ref. [2] (labelled "MV") as compared 
to Kovchegov's model [5] . 



wavelength gluons, where as the MV curve, generated with a hard cutoff, has the 
fewest. 

Finally, we display the effect of the non-Abelian terms on the gluon number 
density. This is most easily seen on a plot of q 2 times the number density (Fig. 6). 3 
From Eq. (5), it is clear that the magnitude of N c \ governs the importance of these 
contributions. In Fig. 6 we have taken N c x to be 0, 50a -2 , and 100a -2 . What we 
see in this plot is a transfer of gluons from low values of q to higher ones. In fact, 
it is possible to show that the area under this plot is conserved [3], and that this 
shifing of gluons from one energy to another is the only effect of the non-Abelian 
contributions. 

In conclusion, we see that imposing a color-neutrality condition on the nucleons 
eliminates the divergent infrared behavior of the two-point vector potential correla- 
tion function in the MV model. Because we have obtained a well-defined expression 
for the gluon number density, we are able to perform a quantitative investigation 
of the features of the MV model. For example, we have shown that the MV model 
predicts a gluon number density which is proportional to l/q + to all orders in the 



3 ) The fact that this curve approaches a constant for large q 2 demonstrates the 1/q 2 behavior 
described in the text. 
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FIGURE 5. Plot of the gluon number den- 
sity at fixed q + . The nucleon size parame- 
ters a and Aqcd have been fixed so that these 



functions match for q 



oo. We have cho- 



sen N c x — 50a 2 . Plotted are the results of 
Ref. [2] (labelled "MV") as well as results us- 
ing Kovchegov's model [5] . 
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FIGURE 6. Plot of q 2 times the gluon num- 
ber density at fixed q + for various values of the 
color charge density within the uniform version 
of Kovchegov's model [5]. As \ increases, the 
importance of the non-Abelian contributions 
increases. Conversely, \ ~~ * corresponds to 
the Abelian limit. 



charge density. We are able to compute the gluon structure function in this frame- 
work, including its absolute normalization. The details of these and other related 
results may be found in Ref. [3]. 
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